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MEASUREMENTS  IN  REGIONS  OF  SHOCK  WAVE/TURBULENT 
BOUNDARY  LAYER  INTERACTION  FROM  MACH  4  TO  10  FOR  OPEN  AND 
“BLIND”  CODE  EVALUATION/VALIDATION 


ABSTRACT 

Experimental  studies  have  been  conducted  to  examine  the  characteristics  of  regions  of  shock 
wave/turbulent  boundary  layer  interaction  over  cone/flare  and  hollow  cylinder/flare  configurations  in  high 
Reynolds  number  flows  at  Mach  numbers  between  5  and  10  in  cold  flows  and  in  flows  with  duplicated 
flight  velocities.  Detailed  surface  pressure  and  heat  transfer  measurements  have  been  made  through  the 
separated  interaction  regions  as  well  as  upstream  and  downstream  of  the  interactions  for  a  range  of 
Reynolds  numbers  under  cold  wall  conditions  for  total  enthalpies  duplicating  flight  The  large  scale  of  the 
models  used  in  these  experiments  enabled  us  to  obtain  measurements  in  fully  turbulent  flows  with  the 
length  of  turbulent  flow  up  to  1 ,000  boundary  layer  thicknesses  downstream  of  the  beginning  of  untripped 
transition.  The  surface  measurements  obtained  in  these  studies  together  with  Schlieren  and  interferometry 
measurements  of  the  regions  of  shock  wave/boundary  layer  interaction  have  been  assembled  to  provide 
data  sets  for  “blind”  code  validation  studies.  We  have  also  included  detailed  information  on  the  results  of 
experimental  studies  of  shock  wave/boundary  layer  interaction  conducted  earlier  which  can  also  be  used  to 
evaluate  the  models  of  turbulence  employed  in  CFD  codes. 

INTRODUCTION 

The  accurate  prediction  of  the  scale  and  aerothermal  loads  generated  in  separated  regions  of  shock 
wave/turbulent  boundary  layer  interaction  remains  one  of  the  most  formidable  computational  tasks, 
principally  because  of  the  difficulty  in  modeling  the  nonequilibrium  turbulence  and  shock  interaction  flow 
phenomena  in  the  separation  and  reattachment  regions.  To  predict  these  flows  with  RANS  techniques,  the 
nonequilibrium  nature  of  the  turbulent  processes  in  the  separation  and  reattachment  regions  requires  models 
that  are  significantly  more  sophisticated  than  those  required  to  describe  attached  flows  in  mild  pressure 
gradients.  When  more  sophisticated  numerical  techniques  are  employed  to  describe  these  flows,  such  as 
LES,  again  the  “Achilles  heel”  is  the  modeling  of  the  wall  layer  and  the  unsteady  shock-  turbulence 
interaction  phenomena,  which  are  key  to  describing  the  mechanisms  of  flow  separation  and  reattachment. 
In  order  to  advance  the  development  of  both  RANS  and  LES  methods,  it  is  of  major  importance  to  provide 
detailed  measurements  in  fully  turbulent  separated  regions  over  well-defined  configurations  in  both  “true 
temperature”  and  “cold”  supersonic  and  hypersonic  flows  with  boundary  and  flow  conditions  to  provide  the 
basis  for  selecting  and  evaluating  the  models  of  turbulence.  Currently  there  are  virtually  no  measurements 
under  cold-wall  conditions  where  the  total  enthalpy  of  the  flow  matches  the  flight  Mach  number  that  are 
required  to  resolve  key  questions  associated  with  compressibility  and  real  gas  effects  on  the  size  and 
structure  of  separated  turbulent  interaction  regions.  The  experiments  conducted  in  the  current  studies  which 
were  made  at  total  enthalpies  which  corresponded  to  flight  velocity  at  the  test  Mach  number  between  5 
and  9  were  designed  to  remedy  this  situation. 

Most  of  the  experiments  conducted  over  the  past  50  years  to  provide  data  on  both  wedge-  and  shock- 
induced  turbulent  separated  regions  in  supersonic  and  hypersonic  flows  suffer  from  deficiencies  associated 
with  the  transitional  nature  of  the  flow  upstream  of  the  interaction  and  well-defined  boundary  conditions 
associated  with  the  experimental  configuration  or  extraneous  disturbances  introduced  into  the  tunnel  wall 
boundary  layer  upstream  of  the  interaction.  There  continues  to  be  a  debate  as  to  whether  the  gross  unsteady 
nature  of  the  separated  region  observed  in  many  experiments  results  from  transition  or  tripping  on  models 
supported  in  the  tunnel  or  from  upstream  effects  on  the  tunnel  walls.  While  many  experiments  with 
separated  interaction  regions  have  been  conducted  with  two-dimensional  configurations,  the  lack  of  well- 
defined  boundary  conditions  has  always  provided  a  major  stumbling  block  to  the  numerical  simulation  of 
these  flows.  This  is  true  even  with  the  large  flat-plate  wedge-  and  flat-plate  shock-generated  configurations 
that  have  been  employed  to  study  laminar  and  turbulent  shock  boundary  layer  interaction  over  a  large  range 
of  test  conditions  at  CUBRC.  In  this  experimental  program  we  have  employed  axisymmetric 
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configurations  such  as  hollow  cylinder/flare  and  double  cone  models  to  provide  measurements  on 
configurations  with  well-defined  boundary  conditions.  In  these  experiments,  the  regions  of  boundary  layer 
transition  are  well  defined  and  are  well  upstream  of  the  interaction  regions.  Also,  the  conditions  in  the 
freestream  are  verified  both  by  surveys  and  computations  from  the  reservoir  conditions. 

The  major  objective  of  the  current  program  is  to  provide  detailed  surface  measurements  in  regions  of 
shock  wave/turbulent  boundary  layer  interaction  to  produce  datasets,  for  “blind”  code  validation  activity 
similar  to  that  we  conducted  earlier  for  laminar  flows.  We  have  also  selected  from  our  earlier  studies  a 
number  of  test  cases  for  which  we  have  presented  the  full  datasets  to  provide  “open  “  measurements  of  the 
distributions  of  heat  transfer  and  pressure  for  the  initial  evaluation  of  the  codes. 

REVIEW  OF  RECENT  COMPARISONS  BETWEEN  PREDICTIONS  AND 
EXPERIMENTS  IN  HYPERSONIC  TURBULENT  INTERACTIONS  (  OPEN 
TEST  CASES) 

Our  recent  experimental  programs  have  employed  large  axisymmetric  models  to  provide  separated 
regions  of  shock  wave/turbulent  boundary  layer  interactions  with  well-defined  model  and  flow  boundary 
conditions,  and  we  have  performed  the  testing  at  high  Reynolds  numbers  to  place  boundary  layer  transition 
at  least  100  boundary  layer  thicknesses  ahead  of  the  shock  interaction  region.  The  large  cone/flare  model 
shown  in  Figure  1(a)  was  constructed  not  only  to  generate  high  Reynolds  number  in  the  interaction  region 
but  also  to  provide  large  boundary  layer  thicknesses  which  could  be  surveyed  with  probes  which  could 
resolve  details  of  the  wall  region  of  the  turbulent  boundary  layer.  Although  the  full-scale  HIFiRE  1  model 
shown  in  Figure  1(b)  was  restricted  in  size  because  of  the  dimensional  requirements  associated  with  the 
launch  system,  measurements  were  obtained  with  well-defined  turbulent  boundary  layers  and,  for  the  cases 
we  have  selected  for  code  validation,  we  achieved  transition  well  upstream  on  the  fore-cone  without 
boundary  layer  trips.  The  length  of  the  flare  again  was  selected  based  on  launch  vehicle  requirements  and 
as  shown  later  was  just  sufficient  to  provide  a  good  downstream  boundary  condition. 


(a)  Large  Cone/Flare  Model  (b)  HIFiRE  1  Model 

Figure  1  Cone  Flare  and  Cone  Cylinder  Flare  Models  to  Examine  Turbulent  Separated  Flows  over 
Cone  Flare  Junctions 

Although  we  obtained  measurements  on  the  large  cone/flare  model,  in  an  experimental  program  before 
we  tested  the  HIFiRE  1  configuration,  at  that  time  we  did  not  have  the  resources  to  compare  data  with 
RANS  calculations  incorporating  contemporary  turbulence  modeling  as  was  done  during  and  after  our  tests 
with  the  HIFiRE  configuration.  The  initial  flare  configuration  for  HIFiRE  1  was  a  36°  flare  and  it  was  clear 
from  our  initial  tests  duplicating  flight  conditions  that,  although  pretest  predictions  indicated  a  small 
separated  region,  our  measurements  indicated  an  extensive  separated  flow  spreading  from  the  downstream 
edge  of  the  flare  over  the  cylindrical  section  of  the  model  was  developed  which  would  not  provide  a 
satisfactory  code  validation  dataset.  Based  on  subsequent  measurements,  we  selected  a  configuration  with  a 
33°  flare  to  provide  a  well-defined  separated  flow  at  the  cylinder/flare  junction  with  sufficient  distance 
downstream  to  provide  well-defined  downstream  conditions.  For  the  conditions  at  which  these  experiments 
were  conducted,  code  predictions  of  the  flow  upstream  and  in  the  interaction  region  were  obtained  with  the 
Shear  Stress  Transport  (SST)1  Model  and  the  Spalart-Allmaras  (SA)  Model.  Predictions  of  the  turbulent 
flow  upstream  of  the  interaction  regions  (as  shown  in  Figure  2)  surprisingly  gave  very  different  results  and 
clearly  issues  associated  with  the  discrepancies  between  theory  and  experiment  should  be  solved  before 
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embarking  on  comparisons  in  the  interaction  region.  This  discrepancy  between  theory  and  measurement  for 
the  turbulent  heating  levels  over  a  cone  was  observed  earlier  in  experiments  performed  at  Mach  10  flight 
conditions.  Shown  in  Figure  3,  while  the  heat  transfer  measurements  in  laminar  flow  were  in  excellent 
agreement  with  prediction,  measured  turbulent  heating  was  significantly  less  than  predicted  with  the  SA, 
SST  and  Baldwin-Lomax  turbulence  models.  The  calculations  of  the  size  of  the  separated  interaction  region 
(shown  in  Figure  4)  demonstrate  that  the  SST  model  significantly  overpredicts  the  size  of  the  separated 
region,  while  the  SA  model  predicting  an  attached  flow  in  the  comer. 


Figure  2  Comparisons  between  Measurements  Upstream  of  the  Interaction  Region  and 
Predictions  for  a  Range  of  Turbulence  Models 


fr-tt,  T  cone  installed  in  LENS  I  facility 


Run  20 


Figure  3  Comparisons  between  Measurement  Sets  of  Heat  Transfer  at  Mach  10  Flight  Conditions 
and  RANS  Prediction  Methods 


EXPERIMENT  (Run  14) 


Shear  Stress  Transport 
(SST)  Model 
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Figure  4  Comparison  between  Measured  Size  of  Separated  Region  on  HIFiRE  1  Model  and 
Original  Predictions  with  DPLR  using  SST  and  SA  Turbulence  Models 

An  empirical  modification  to  these  methods  suggested  by  Wilcox2  was  employed  to  limit  the 
generation  of  shear  stress  by  introducing  an  empirical  coefficient  CLIM  in  the  shear  stress  production 
equation.  Figure  5  shows  the  results  of  employing  different  values  for  this  coefficient.  For  this  set  of 
experimental  conditions,  a  value  of  CLIM  =  0.9  provides  the  most  accurate  method  of  predicting  the 
separated  flow.  For  this  value  RANS  calculations  are  in  reasonable  agreement  with  the  pressure  and  heat 
transfer  measurements  obtained  in  the  interaction  regions  over  the  cylinder/flare  as  shown  in  this  figure. 
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Figure  5  Comparison  between  DPLR/RANS  Solutions  and  Heat  Transfer  and  Pressure 
Measurements  at  the  Cylinder  Flare  Junction  on  the  HIFiRE  1  Configuration 


Comparisons  between  Computations  and  Measurements  obtained  Earlier  in  Separated  Flow  over 
Flat  Plate/Wedge  and  Flat  Plate/  Shock  Generator  Models 

Earlier  extensive  sets  of  measurements  were  made  with  flat  plate/wedge  and  flat  plate/shock  generator 
configurations  at  CUBRC  over  the  past  30  years.  We  have  selected  two  test  cases  from  this  dataset  which 
we  believe  present  the  most  useful  from  the  viewpoint  of  code  validation.  Two  of  the  models  used  in  these 
studies  are  shown  in  Figure  6. 


(a)  Wedge  Interaction  Model  (b)  Shock  Interaction  Model 

Figure  6  Two-dimensional  Models  to  Examine  Shock  Wave/Turbulent  Boundary  Layer  Interaction 
in  Compression  Corners  and  Externally  Developed  Shocks  Incident  on  the  Turbulent 
Boundary  Layer 

Regions  of  flow  separation  either  induced  at  compression  comers  or  by  incident  shocks  in  hypersonic 
flows  are  complex  in  nature,  principally  because  flow  separation  takes  place  initially  at  the  base  of  the 
turbulent  boundary  layer  and  the  recirculation  region  is  often  embedded  in  the  incoming  boundary  layer  as 
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illustrated  in  the  sequence  of  photographs  for  wedge-  and  shock-induced  separated  regions  shown  in 
Figure  7.  Of  significant  interest  here  is  that  the  wedge-shock  extends  into  the  sublayer  of  the  turbulent 
flows,  initially  creating  a  tiny  embedded  recirculating  region  at  the  flat  plate  wedge  junction  or  at  the  base 
of  the  incident  shock.  In  both  cases,  separation  shock  traverses  and  interacts  with  the  turbulence  in  the 
incoming  boundary  layer.  Only  when  turning  angles  of  over  36°  are  imposed  on  the  flow  are  well- separated 
flows  generated  which  extend  the  full  width  of  the  boundary  layer.  Modeling  the  turbulence  generated 
through  the  interaction  region  of  these  flows  is  clearly  a  formidable  task  and  selecting  the  models  which 
most  accurately  represent  these  phenomena  will  require  detailed  flowfield  as  well  as  surface  measurements. 


(d)  WEDGE  ANGLE  -  36° 


U)  SHOCK  GENERATOR  ANGLE  -  125° 


(b)  SHOCK  GENERATOR  ANGLE  -  15° 


(c)  SHOCK  GENERATOR  ANGLE  -  17.5° 


Figure  7  Schlieren  Photographs  of  Turbulent  Separated  Flows  at  Wedge  Compression  Corners  and 
Induced  by  Incident  Shock/Turbulent  Boundary  Layer  Interaction 

When  we  employed  the  DPLR  method  with  a  turbulence  model  which  best  fit  the  measurements  over 
the  HIFiRE  configuration  to  describe  flow  for  these  two  configurations,  we  significantly  overpredicted  the 
size  of  the  interaction  region.  Figure  8  shows  the  comparisons  for  the  flat  plate/3 6°  wedge  configuration, 
and  Figure  9  show  similar  comparisons  for  the  separated  flows  introduced  by  incident  shock.  In  both  cases, 
the  computations  capture  with  good  accuracy  the  pressure  and  heating  levels  both  upstream  and 
downstream  of  the  separated  interaction  regions.  Clearly,  the  turbulence  models  employed  here  are  totally 
inadequate  in  describing  the  length  and  properties  of  the  interaction  regions  whether  the  pressure  rise  is 
induced  internally  as  in  the  case  of  the  compression  comer  or  externally  by  an  incident  shock.  It  should  be 
noted  here  that  computations  for  the  incident  shock  cases  are  made  significantly  more  complicated  because 
of  the  requirement  to  compute  the  flow  over  the  shock  generator  including  the  position  of  transition  all  of 
which  can  influence  the  position  of  the  incident  shock.  This  added  complexity  has  coupled  with  those 
associated  with  lateral  boundary  conditions  make  the  axisymmetric  cone/flare  and  hollow  cylinder  flare 
models  far  more  suitable  for  code  validation 
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Figure  8  Comparisons  between  DPLR/RANS  Predictions  and  Measurements  of  Pressure  and 
Heat  Transfer  in  Separated  Flows  over  a  Flat  Plate  Wedge  Configuration 
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Figure  9  Comparison  between  DPLR  Navier/RANS  Predictions  and  Measurements  of  Pressure 
and  Heat  Transfer  in  Separated  Flow  Induced  by  an  Externally  Generated  Incident 
Shock  on  a  Flat  Plate 


Next,  we  show  comparisons  between  the  DPLR  calculations,  employing  the  turbulence  models  that 
were  used  “successfully”  to  predict  the  HIFiRE  test  cases,  with  measurements  on  a  9-ft  long  cone/flare 
configuration  constructed  many  years  ago  over  which  the  boundary  layer  is  fully  turbulent  over  90  percent 
of  the  model  length.  Comparisons  between  the  RANS  predictions  with  the  SA  model  and  Wilcox  stress 
limiter  and  experiment,  shown  in  Figure  10,  indicate  that  the  separated  region  was  significantly 
overpredicted,  although  the  levels  of  pressure  and  heating  upstream  and  downstream  of  the  interaction 
region  are  relatively  accurately  predicted. 
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distance  from  nosetip,  inches 


Figure  10  Comparison  between  DPLR  Calculations  and  Measurements  of  the  Pressure  and  Heat 

Transfer  Distribution  at  the  Cone/Flare  Junction  on  the  Large  Cone  Model 


Diagrams  of  the  models  and  instrumentation  and  tabulations  of  the  measurements  made  for  these  test 
cases  are  presented  in  Appendices  A.l-A-4 


MEASUREMENTS  FOR  “BLIND”  CODE  VALIDATION  STUDIES  IN  REGIONS  OF 
SHOCK  WAVE/TURBULENT  BOUNDARY  LAYER  INTERACTION  ON  NEW 
CONE/FLARE  AND  HOLLOW  CYLINDER/FLARE  CONFIGURATION  AT  MACH 
NUMBERS  FROM  6  TO  9  AT  FLIGHT  MATCHED  ENTHALPIES 

INTRODUCTION 

We  have  performed  experiments  with  a  large  double  cone  and  hollow  cylinder/flare  model  to  examine 
regions  of  shock  wave/turbulent  boundary  layer  interaction  for  a  range  of  high  Reynolds  number  flow 
conditions  from  Mach  5  to  8.  We  will  use  these  as  a  basis  of  “blind”  validation  exercise  where  we  will 
provide  only  details  of  the  model  geometries  and  the  test  conditions  at  which  these  models  were  conducted. 
The  models  were  instrumented  to  provide  highly  detailed  surface  pressure  and  heat  transfer  measurements 
both  upstream  and  downstream,  as  well  as  in  the  interaction  regions  at  the  base  of  the  flares.  A  key  feature 
of  these  boundary  layer  studies  is  that  the  models  employed  in  this  program  are  very  large  so  that  transition 
occurs  well  upstream  of  the  interaction  region,  and  thus  the  influence  of  transition  on  the  unsteady 
characteristics  of  the  interaction  region  and  the  requirement  to  accurately  predict  the  characteristics  of  the 
transition  zone  were  minimized.  However,  we  have  also  obtained  detailed  measurements  in  the  transition 
regions  to  evaluate  the  performance  of  the  prediction  methods  in  these  regions. 

EXPERIMENTAL  FACILITIES 

The  experimental  studies  reported  here  have  been  conducted  in  the  CUBRC  reflected  shock  tunnels. 
A  summary  of  the  performance  of  CUBRC’ s  three  major  ground  test  facilities,  the  LENS  I  and  II  shock 
tunnels,  and  LENS  XX  expansion  tunnel  are  listed  in  Figure  11.  Velocity/altitude  performance  of  the  three 
LENS  facilities  is  presented  in  Figure  12  and  the  Mach  Number/Reynolds  Number  Capabilities  for  LENS  I 
and  II  are  shown  in  Figure  13.  The  LENS  II  tunnel,  which  has  been  employed  in  most  of  the  studies  of 
shock  wave/turbulent  boundary  layer  interaction,  is  operated  principally  to  examine  flows  at  Mach  numbers 
between  3.5  and  10,  and  it  is  capable  of  fully  duplicating  test  conditions  on  vehicles  up  to  30  ft  in  length.  In 
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this  tunnel,  we  also  employ  electrically  heated  models  to  bring  the  surface  temperature  of  the  vehicle  to 
flight  values.  The  LENS  I  facility  overlaps  with  LENS  II  between  Mach  8  and  10  and  can  be  used  to 
duplicate  velocities  up  to  15,000  ft/s  and  Mach  numbers  up  to  22. 


LENS  I 

Velocity  Range  (ft/sec) 
3,000-15,000 
Altitude  (k  ft) 
25-300 
Mach  Numbers 
7-16  (true  duplication) 
16-24 (Re-M  only) 
Reynolds  Numbers (1 /ft) 
104-108 
Test  Time  (ms) 
up  to  25  ms 
Nozzles 

Mach  7-9  (44" Exit) 
Mach  9  -  16  (44"  Exit) 
Mach  15-24  (72"  Exit) 


LENS  II 

Velocity  Range  (ft/sec) 
2,500-10,000 
Altitude  (k  ft) 

SL -  200 


LENS XX 

Velocity  Range  (ft/sec) 
8,000-40,000 
Altitude  (k  ft) 

SL -  250 


Mach  Numbers 


Mach  Numbers 


3-10 

Reynolds  Numbers  (1  /ft) 
105-109 


8-25 

R  ey  n  olds  N  u  mbers  ( 1 /ft) 
103-109 


Test  Time  (ms) 
up  to  100  ms 
Nozzles 


Test  Time 
up  to  4  ms 
Nozzles 


Mach  3 -5  (42"  Exit) 
Mach 5 -7  (60” Exit) 
Mach  7  - 10  (72"  Exit) 


Contoured  for  Mach 
Numbers  from  8  to  25 


Figure  11  CUBRC  Test  Facilities  to  Investigate  Shock  Wave/Turbulent  Boundary  Layer 
Interaction  Phenomena  in  Supersonic  and  Hypervelocity  Flows 
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Figure  12  LENS  Facilities  Altitude/Velocity  Map 
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— Performance  Profile  of  LENS  I  — Performance  Profile  of  LENS  II  — Performance  Profile  of  48*  Tunnel 

a  LENS  I  D-Noule  Calibrations  o  LENS  I C  Nozzle  Calibrations  *  LENS  II  M4  Nozzle  Calibrations 

o  LENS  II  M6  No?7le  Calibrations  o  48*  Tunnel  Calibrations  — 3  5%  Scale  Space  Shuttle  (OTS) 


Figure  13  Mach  Number/Reynolds  Number  Envelope  for  LENS  Tunnels 

The  experimental  studies  conducted  in  the  LENS  II  tunnel  were  conducted  at  freestream  Mach 
numbers  between  4  and  9  and  Reynolds  numbers  up  to  0.5  billion.  They  were  also  conducted  for  both 
“cold  flow”  and  at  flight  enthalpies  at  each  Mach  number.  The  large  scales  of  the  models  which  are 
employed  in  these  studies,  coupled  with  the  ability  to  generate  Reynolds  numbers  up  to  80  x  1 06  per  foot, 
enables  us  to  obtain  transition  within  several  inches  of  the  leading  edge  of  the  models,  thereby  eliminating 
any  questions  as  to  the  unsteady  effects  of  transition  on  the  mean  and  fluctuating  characteristics  of  the 
interaction  regions.  In  turn,  the  specific  models  employed  to  describe  the  transition  regions  in  the  CFD 
codes  also  must  have  little  effect  on  the  predictions  of  the  characteristics  of  the  shock  interaction  regions. 
During  these  studies,  we  have  employed  both  intrusive  and  nonintrusive  measurements  to  define  the  mean 
and  fluctuating  characteristics  of  the  freestream.  Using  this  information,  we  are  evaluating  the  methods  to 
predict  transition  onset  (i.e.,  STABL4,  the  CRAFT  ETM5,  and  semi-empirical  design  methods)  in  our 
ground  test  environment  with  the  objective  of  developing  rational  methods  to  link  ground  test 
measurements  to  flight. 

Evaluation  of  Test  Conditions 

The  facility  test  conditions  are  established  and  validated  by  a  combination  of  measurements  in  the 
reservoir  region  and  test  section  of  the  tunnel.  The  stagnation  pressure  and  enthalpy  are  obtained  directly 
from  pressure  measurements  behind  the  reflected  shock  and  measurements  of  the  incident  shock  Mach 
number.  Conditions  in  the  freestream  are  obtained  from  survey  rakes  containing  pitot  pressure  probes, 
stagnation  heat  transfer  gauges  on  hemispherical  cylinders,  total  temperature  measurements  with  vented 
thermocouple  probes  (where  applicable),  and  cone  pressure  measurements  to  establish,  in  conjunction  with 
the  aforementioned  measurements,  the  static  pressure  in  the  flow.  Facility  calibration  is  always  supported 
by  comparisons  to  numerical  predictions.  Calculations  of  the  nozzle  flowfield  at  each  select  test  conditions 
have  been  performed  using  a  specialized  form  of  the  DPFR  code6.  It  employs  the  Spalart-Allmaras  1- 
equation  turbulence  model  (Spalart  and  Allmaras7)  with  the  Catris  and  Aupoix8  compressibility  correction. 
Non-ideal  gas  effects  are  included  in  the  calculation  because  of  the  large  pressures  in  the  reservoir.9  From 
these  measurements  and  calculations,  we  can  determine  the  accuracy  of  the  freestream  dynamic  pressure 
and  the  stagnation  point  enthalpy  to  ±  5%  and  the  Mach  number  to  ±1.5%. 

MODELS  AND  INSTRUMENTATION 

Data  sets  for  the  “blind”  validation  studies  were  obtained  with  high-frequency  surface  pressure  and 
heat  transfer  instrumentation  installed  in  the  new  large  cone/flare  model  and  the  large  hollow  cylinder/flare 
model  shown  installed  in  the  72-inch  nozzle  in  the  schematic  diagrams  in  Figure  14.  The  double  cone 
model  was  tested  first  in  the  FENS  II  tunnel  in  our  Mach  4.5-to-6  nozzle  with  a  60-inch  exit  plane  diameter 
as  shown  in  Figure  15.  The  dimensions  of  this  model  are  presented  in  Figure  16.  The  model  is  instrumented 
along  its  length  with  high-frequency  thin-film  heat  transfer  and  piezoresistive  pressure  sensors,  and  as 
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shown  in  Figure  17,  there  is  a  high  density  of  instrumentation  at  the  cone/flare  junction.  In  addition  to 
Schlieren  photography,  flow  field  measurements  were  made  with  differential  interferometry. 


Figure  14  Large  Hollow  Cylinder/  Flare  and  Cone/Flare  Model  installed  in  the  72inch  Mach  5  to  10 

Contoured  Nozzle 


Figure  15  New  Large  Cone/Flare  Model  Constructed  and  Employed  in  the  LENS  II  Tunnel  at 
Mach  6  to  8  to  Generate  New  Datasets  for  Blind  Code  Validation  Studies 


Figure  16  Model  Dimensions  and  Instrumentation  Layout  of  the  New  Large  Cone/Flare  Model 
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Figure  17  Instrumentation  at  Cone/Flare  Junction  on  Large  Cone  Model  and  Schlieren 

Photograph  of  Transitional  Flow  over  the  Cone/Flare  Junction  (Not  a  new  Test  Case) 

The  dimensions  of  the  new  hollow  cylinder/flare  model  are  shown  in  Figure  18.  This  model  is  also 
highly  instrumented  at  the  cylinder/flare  junction  to  provide  detailed  surface  measurements  on  the  size  and 
turbulence  properties  throughout  the  interaction  region  as  well  as  the  attached  turbulent  flows  both  ahead 
and  downstream  of  the  separated  interaction  region  This  model  is  shown  installed  in  the  exit  plane  of  the 
new  Mach  7  to  10  nozzle  in  Figure  19.  A  photograph  of  the  large  hollow  cylinder/flare  model  is  shown  in 
Figure  20  prior  to  installation  in  the  LENS  II  tunnel. 


109,02 


Figure  18  Model  Dimensions  and  Instrumentation  Layout  of  the  Large  Hollow  Cylinder/Flare 
Model 
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Figure  19  Large  Hollow  Cylinder  Flare  Installed  in  the  LENS  II  72-inch  Contoured  Mach  5  to 

10  Nozzle 


Figure  20  Photograph  of  Large  Hollow  Cylinder  Flare  Model 


NEW  DATASETS  FOR  CODE  VALIDATION 

Datasets  have  been  collected  from  measurements  made  with  both  the  large  cone/flare  model  and 
hollow  cylinder/flare  model  from  tests  conducted  in  the  LENS  facility  at  Mach  numbers  from  6,  7  and  8  at 
a  number  of  high  Reynolds  number  test  conditions.  Table  I  shows  a  typical  set  of  basic  freestream 
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conditions  employed  in  the  Mach  6  studies.  Here  it  can  be  seen  that  measurements  were  obtained  for  both 
cold  flow  conditions  and  for  true  temperature  flight  conditions.  Similar  sets  of  measurements  are  being 
obtained  on  both  the  cone/flare  and  hollow  cylinder/flare  models  at  Mach  numbers  of  7  and  8.  Again,  at 
cold  flow  conditions  and  for  conditions  where  we  matched  total  enthalpy  at  the  specific  Mach  number 
tested. 


Run 

Number 

Mach 

Number 

Velocity 

Temperature 

Pressure 

Density 

Reynolds 
Number 
(Cone  Length) 

Tw/To 

- 

(ft/sec) 

{deg  R) 

(psia) 

(slugs/ftA3) 

(ReL) 

- 

1 

6.04 

3138.90 

112.16 

0.25 

1.85E-04 

5. 15 E +07 

0.57 

4 

5.87 

3249.50 

127.04 

0.12 

7.66E-05 

1.92E+07 

0.53 

5 

6.00 

5467.70 

344.82 

0.54 

1.27E-04 

2.05E+07 

0.20 

9 

7.15 

3485.03 

98.89 

0.05 

4.00E-05 

1.44E+07 

0.48 

10 

7.01 

6277.46 

333.80 

0.34 

8. 67  E -05 

1.62E+07 

0.15 

11 

7.13 

3440.73 

96.93 

0.32 

2. 73  E -04 

9.89E+07 

0.49 

12 

7.00 

6483.02 

357.04 

0.17 

4.05  E -05 

7.35E+06 

0.14 

14 

7.14 

3380.71 

93.32 

0.16 

1.42E-04 

5.28E+07 

0.52 

17 

7.16 

3381.55 

92.85 

0.02 

2. 10 E -05 

7.86E+06 

0.62 

18 

7.11 

3457.93 

98.46 

0.32 

2.76E-04 

9.87E+07 

0.49 

Table  1  Freestream  Condition  for  Existing  Large  Cone  Flare  Validation  Cases  at  Mach 

Numbers  of  6  and  7 


Shown  in  Figure  21  are  heat  transfer  measurements  made  during  the  experimental  program  of  the 
heating  rate  along  the  cone  demonstrating  that  transition  occurs  well  ahead  of  the  interaction  region 
resulting  in  a  well-defined  turbulent  boundary  layer  upstream  of  the  cone/flare  junction. 


7°  Cone  Flare  Heat  Transfer  Distribution 


0  20  40  60  80  100 


Axial  location  (in) 


Figure  21  Heat  Transfer  Measurements  along  the  Cone  Section  of  the  Cone/Flare  Model 
demonstrating  the  Fully  Turbulent  Nature  of  the  Flow  ahead  of  the  Cone/Flare 
Junction 

Ideally,  the  computations  to  describe  the  shock  interaction  region  should  be  initiated  at  the  leading 
edge  of  both  the  hollow  cylinder/flare  and  double  cone  model.  The  calculations  should  describe  accurately 
the  laminar  boundary  layer  downstream  of  the  leading  edge,  the  transition  region  from  laminar  to  turbulent 
flow,  and  the  region  of  shock  wave/turbulent  boundary  layer  interaction.  An  important  boundary  condition 
for  these  calculations  is  the  position  of  boundary  layer  transition.  So  to  specify  this  point  for  each  of  the  test 
cases,  we  have  shown  the  laminar  heating  distribution  on  each  model  up  to  the  point  where  the  heat  transfer 
rate  indicates  boundary  layer  transition  has  occurred.  Shown  in  Figure  22  are  typical  measurements  taken 
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on  the  forecone  of  the  cone/flare  model  defining  the  laminar  heating  levels  and  the  position  at  which  the 
heating  levels  show  that  boundary  layer  transition  is  occurring. 


Run 

Number 

Laminar  Heating 
Value  at  Transition 

Location 

Location  of 

Start  of 

Transition 

Location  of 
Completion  of 
Transition 

(Btu/ftA2-sec) 

(inches) 

(inches) 

1.00 

0.71 

6.85 

13.89 

4.00 

0.28 

15.89 

25.26 

5.00 

6.51 

9.95 

16.89 

9.00 

0.27 

18.89 

28.84 

10.00 

6.42 

17.89 

28.84 

11.00 

Turbulent  Forward  of  First  Sensor 

12.00 

3.34 

27.83 

45.82 

14.00 

7.18 

15.89 

28.84 

17.00 

3.42 

48.52 

74.33 

18.00 

Turbulent  Forward  of  First  Sensor 

Figure  22  Laminar  Heating  Levels  and  Transition  Locations  for  Double  Cone  Studies  with  Turbulent 
Heating 


As  discussed  earlier,  all  the  measurements  for  the  “blind”  code  validation  studies  are  being  obtained  at 
freestream  conditions  where  boundary  layer  transition  occurs  well  upstream  of  the  interaction  regions  at 
conditions  similar  to  those  where  the  heat  transfer  measurements  shown  in  Figure  21  were  made.  However, 
in  companion  studies  where  we  were  specifically  looking  at  the  characteristics  of  the  transition  process, 
measurements  were  made  at  lower  Reynolds  numbers  and  for  these  conditions  the  separated  flow  over  the 
cone/flare  junction  was  transitional  in  nature.  A  Schlieren  photograph  of  this  transitional  flow  is  shown  in 
Figure  23  and  clearly  this  is  not  one  of  the  test  cases  for  the  “blind”  validation  studies;  however,  we  have 
shown  the  measured  pressure  distribution  through  this  region  in  Figure  24  as  an  example  of  the  detailed 
measurements  made  through  the  interaction  region  in  the  subsequent  measurements  at  significantly  larger 
Reynolds  numbers. 
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Figure  23  Transitional  Flow  over  Cone/Flare  Model  at  Much  Lower  Reynolds  Numbers  than 

used  for  the  Fully  Turbulent  Interaction  Studies 
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Figure  24  Distribution  of  Pressure  for  Transitional  Interaction 


SUMMARY  AND  CONCLUSIONS 

Experimental  studies  have  been  conducted  to  examine  the  characteristics  of  regions  of  shock 
wave/turbulent  boundary  layer  interaction  over  cone/flare  and  hollow  cylinder/flare  configurations  in  high 
Reynolds  number  flows  at  Mach  numbers  from  5  to  10  in  cold  flows  and  at  duplicated  flight  velocities. 
Detailed  surface  pressure  and  heat  transfer  measurements  have  been  made  through  the  separated  interaction 
regions  as  well  as  upstream  and  downstream  of  the  interactions  for  a  range  of  Reynolds  numbers  under  cold 
wall  conditions  for  total  enthalpies  matching  flight,  as  well  as  at  lower  enthalpies  which  matched 
measurements  made  earlier  for  low  temperature  freestream  conditions.  The  large  scale  of  the  models  used 
in  these  experiments  enabled  us  to  obtain  measurements  in  fully  turbulent  flows  with  the  length  of  turbulent 
flow  up  to  1,000  boundary  layer  thicknesses  downstream  of  the  beginning  of  untripped  transition.  The 
surface  measurements  obtained  in  these  studies  together  with  Schlieren  and  interferometry  measurements 
of  the  regions  of  shock  wave/boundary  layer  interaction  have  been  assembled  to  provide  data  sets  for 
“blind”  code  validation  studies.  We  have  also  included  detailed  information  on  the  results  of  experimental 
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studies  of  shock  wave/boundary  layer  interaction  conducted  earlier  which  can  also  be  used  to  evaluate  the 
models  of  turbulence  employed  in  CFD  codes. 
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Appendix  A, 


A.1  Flat  Plate/  Wedge  Data  and  Conditions  for  Run  54 


H  mi 

_ j 

r 

*1 

p,  sl/ft3 

T,  R 

U,  ft/s 

M 

1.600x1  O'4 

110 

5,802 

11.3 

X,  in 

qw,  btu/ft2-s 

X,in 

p,  psia 

X,in 

tau,  psia 

32.52 

4.39 

32.53 

0.206 

32.53 

0.00615 

34.02 

4.82 

35.02 

0.221 

36.72 

0.00868 

35.04 

4.5 

36.72 

0.204 

37.21 

0.01 

35.93 

4 

37.31 

0.198 

37.75 

0.00204 

36.13 

4.13 

37.75 

0.327 

38.27 

-0.00299 

36.48 

4.19 

38.27 

0.76 

39.56 

-0.0055 

36.69 

4.1 

38.82 

1.15 

40.48 

0.0114 

37.56 

4 

39.56 

1.08 

42.39 

0.571 

37.76 

3.98 

41.78 

10.2 

43.32 

0.4 

38.08 

9.5 

42.38 

16.5 

45.27 

0.245 

38.61 

8.85 

43.02 

16.6 

47.19 

0.247 

38.82 

10.5 

43.67 

14.3 

39.13 

10.4 

44.31 

14.9 

39.28 

9.81 

44.94 

15.4 

39.42 

10.2 

45.59 

14.9 

39.59 

9.91 

46.22 

14.6 

39.68 

10.9 

47.5  14.5 
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A.2  Shock  Generator/  Flat  Plate  .  Data  and  Conditions  for  Run  49 


^ - 

h - - 50* 


3.000* 

T" 


T 

9.919* 


*1 


p,  sl/ft3 

T,  R 

U,  ft/s 

M 

1.729x1  O'4 

113 

5,950 

11.4 

X,  in 

qw,  btu/ft2-s 

X,in 

p,  psia 

X,in 

tau,  psia 

32.52 

5.16 

32.53 

0.241 

32.53 

0.00662 

34.02 

5.04 

35.02 

0.323 

36.72 

-0.00562 

35.04 

4.25 

36.72 

1.35 

37.21 

-0.00481 

35.93 

10 

37.31 

1.25 

37.75 

-0.05 

36.13 

11.1 

37.75 

5.26 

38.27 

-0.0155 

36.48 

13.2 

38.27 

17.5 

39.56 

0.4 

36.69 

13.6 

38.82 

21.9 

40.48 

0.485 

37.56 

30.4 

39.56 

27.2 

42.39 

0.413 

37.76 

50 

41.78 

26.8 

45.27 

0.181 

38.08 

92.4 

42.38 

27.1 

47.19 

0.127 

38.61 

220 

43.02 

29.4 

38.82 

218 

43.67 

21.5 

39.13 

254 

44.31 

16.4 

39.28 

248 

44.94 

12.2 

39.42 

254 

45.59 

9.33 

39.59 

221 

46.22 

7.72 

39.68 

235 

47.5  4.81 

39.89 

233 

48.79 

3.96 
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A.3  Large  Cone  Flare  Model  Data  and  Conditions  for  Large  Cone/Flare  Run  4 


1  U4 

T 

-  1U  —*\ 

p,  sl/ft3 

T,  R 

U,  ft/s 

M 

6.345xl0'5 

121 

5,920 

11.0 

X,  in 

CH-q/[pU(h0-hw)] 

X,in 

p/qoo 

98.91 

0.000805 

98.910 

0.034617 

101.32 

0.000787 

100.12 

0.037367 

101.61 

0.000798 

101.92 

0.044641 

101.92 

0.000847 

102.51 

0.104292 

102.22 

0.001240 

103.11 

0.167478 

102.51 

0.002023 

103.72 

0.179795 

102.81 

0.002017 

104.02 

0.198549 

103.11 

0.001816 

104.68 

0.168566 

103.42 

0.001830 

106.18 

1.931486 

103.72 

0.001642 

106.68 

1.957518 

104.02 

0.002398 

107.18 

1.639943 

105.18 

0.012382 

107.68 

1.812071 

105.68 

0.023679 

108.68 

1.880974 

106.18 

0.027989 

109.68 

1.997410 

106.68 

0.027113 

110.68 

1.677503 

107.18 

0.024481 

111.68 

1.553167 

107.68 

0.024153 

112.68 

1.590079 
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A.4  HIFiRE-1/2  Cone/Cylinder/Flare  Run  30 


270* 


p,  sl/ft3 

T,  R 

U,  ft/s 

M 

1.3X10-4 

408 

7,120 

7.19 

X,  in 

qw,  btu/ft2-s 

X,in 

p,  psia 

6.848 

18.920 

26.834 

1.553 

9.892 

15.984 

39.342 

1.561 

13.892 

12.151 

48.516 

0.503 

14.892 

12.954 

53.266 

0.485 

15.892 

12.384 

57.516 

0.479 

16.892 

12.357 

60.643 

0.456 

17.892 

18.669 

60.953 

0.494 

18.892 

24.482 

61.108 

0.424 

19.892 

26.909 

61.263 

0.754 

20.892 

36.284 

61.418 

0.780 

21.892 

41.865 

61.573 

0.749 

22.761 

47.132 

61.728 

0.755 

23 


63.943 


48.892 


64.073 

59.491 

64.203 

71.687 

64.333 

85.973 

64.463 

107.734 

64.593 

118.448 

64.723 

139.267 

64.853 

156.220 

64.983 

176.113 

65.113 

205.019 

65.243 

228.316 

65.373 

239.008 

65.503 

266.669 

65.633 

293.006 

65.763 

304.615 

65.893 

273.236 

66.368 

265.986 

66.535 

255.345 

66.703 

259.736 

66.871 

255.090 

67.039 

254.170 

67.206 

247.960 

67.374 

244.740 

67.542 

249.010 

67.710 

236.400 

X  represents  axial  distance  from  the  virtual  sharp  nose  (2.5-mm  nose  radius  for  Run  30) 


